﻿//struct CoutteFlowIC : WorkletInitCondition<CoutteFlowIC>
//{
//  VTKM_EXEC ValueType ComputeValue(const Vec3f& p) const
//  {
//    FloatDefault gam = 1.4, gamm1 = gam - 1, mach = 0.2;
//    FloatDefault u = 1;
//    FloatDefault v = 0;
//    FloatDefault density = 1;
//    FloatDefault pre = 1. / (gam * mach * mach);
//
//    ValueType value;
//    value[0] = density;
//    value[1] = density * u;
//    value[2] = density * v;
//    value[3] = density * 0;
//    value[4] = pre / gamm1 + 0.5 * density * (u * u + v * v);
//
//    return value;
//  }
//};